Three-dimensional Binary Superlattices of Oppositely-charged Colloids 
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We report the equilibrium self-assembly of binary crystals of oppositely-charged colloidal micro- 
spheres at high density. By varying the magnitude of the charge on near equal-sized spheres we 
show that the structure of the binary crystal may be switched between face-centered cubic, cesium 
chloride and sodium chloride. We interpret these transformations in terms of a competition between 
entropic and Coulombic forces. 



The spontaneous crystallization of mixtures of large 
and small colloids continues to fascinate material scien- 
tists since the first observation of colloidal alloys, over 
two decades ago 0. There are several reasons for this 
continuing interest: From a fundamental standpoint, col- 
loids provide arguably the best experimental realization 
of the hard sphere model, whose phase diagram is com- 
pletely dominated by entropy. Consequently the rich va- 
riety of self-assembly phenomena seen in colloidal mix- 
tures provides a fascinating test bed for many-body sta- 
tistical physics. From a practical standpoint, colloidal 
assembly has great technological potential as a fabrica- 
tion technique for three-dimensional photonic band gap 



(PBG) materials 
chemical sensors 



, ultrafast optical switches [4| and 
Most studies of colloidal assembly 
have, to date, been restricted to the simplest situation, 
that of purely repulsive spherical particles where the par- 
ticles bear either no or similar surface charges and so be- 
have as hard spheres. Extensive computer simulations 
0, and experiments have located 0, 0, E El El only 
three equilibrium binary phases for hard spheres, AB, 
AB2 and AB13, iso-structural with the atomic analogues 
NaCl, AIB2 and NaZni3. The small number and the 
essentially close-packed nature of these lattices limit the 
practical usefulness of colloidal assembly techniques, par- 
ticularly in the search for a visible PBG material. More 
complex architectures could, in principle, be achieved by 
introducing weak reversible attractions betwee n sp ecies 
or by fine-tuning the inter-particle forces [3 El EH Un- 
fortunately all attempts to date to direct the assembly of 
attractive binary particles have created highly-disordered 
aggregates rather than the hoped-for crystalline struc- 
tures. 

In this letter, we demonstrate how the introduction of 
a weak reversible attraction between unlike particles sta- 
bilizes non-close packed superlattices which would other- 
wise be entropically unfavorable. We use fluorescent con- 
focal microscopy to follow the self-organization of a non- 



bc 



addressed 



tTo whom correspondence should 
(p.bartlett@bristol.ac.uk) 

^Current address: Soft Condensed Matter Group, Debeye Institute, 
University of Utrecht, The Netherlands. 



aqueous binary suspension of micron-sized spheres with 
similar diameters but carrying opposite charges. The 
non-polar environment used ensures that the strength of 
the attraction is comparable to fc^T and all interactions 
remain reversible. We study the phase behavior as a func- 
tion of the relative strength of the attractive and thermal 
interactions and find that we can tune the translational 
crystal symmetry from face-centered cubic (fee), through 
body-centered cubic (bee) to simple cubic (sc). 

Our system consists of a binary mixture of near- 
identical sized colloidal poly(methyl methacrylate) par- 
ticles [16] suspended in a non-polar density and index- 
matched solvent. Particles are distinguished by con- 
taining one of two different cationic fluorescent dyes. 
The large spheres (component A) have a core radius 
of Ra = 777 nm and contain the orange-red fluores- 
cent dye DilCis (maximum emission at A em = 565 nm). 
The small spheres (component B), with a core radius of 
Rb = 720 nm, are labeled with the green fluorescent dye 
DiOCig (A em = 502 nm). The radius ratio (7 = Rb/Ra) 
is accordingly very close to unity (7 = 0.93 ±0.01). Both 
particles have a narrow distribution of sizes (polydisper- 
sity a of about 4%) and individually form random hexag- 
onal close-packed (rhep) crystals at high densities. The 
volume fraction (f> of the samples was defined relative to 
the hard-sphere freezing transition at 0/ = 0.494. The 
solvent consists of a 78:22 wt % density and near-index 
matching mixture of cycloheptyl bromide (CHPB) and 
cis-decalin (CD). PMMA particles dispersed in this low 
polarity solvent develop a small positive surface charge. 
The charge on each particle is controlled by the polar- 
ity of the dye used as the fluorescent label, with parti- 
cles containing the polar DilCis dye developing a larger 
charge than spheres labeled with the less polar DiOCis 
molecule. The particle charge may be reduced by increas- 
ing the Br - ion concentration in solution jl7j . Generat- 
ing free Br~ by irradiating the suspension with ultravio- 
let light or by including a small length of ferromagnetic 
wire to act as a site for catalytic decomposition of cy- 
cloheptyl bromide, led to a reduction in particle charge 
until at the point of zero charge (PZC) the sign of the 
charge reversed and then became progressively more neg- 
ative with further increases in [Br - ]. Charge inversion 
provides a useful, and to date unexplored, tool to ad- 
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just the inter-particle interactions in a binary suspen- 
sion. By working near the point of zero charge the size 
and sign of the charge on each species may be system- 
atically controlled and the total inter-particle potential 
"tuned" . The point of zero charge is a function of the 
initial particle charge and so is different for DilCis an d 
DiOCis-labeled particles. At low (or high) Br~ concen- 
trations both particles are positive (or negative), all inter- 
actions repulsive and a hard sphere model is appropriate. 
However, between these two limits the two spheres have 
opposite signs and the unlike interactions are attractive. 
We found that, in our system, we could achieve partial 
charge inversion by adding a small length of ferromag- 
netic wire to our suspensions. The charge on the particle 
is a function of the time, t, of contact between suspen- 
sion and metal. At t = both particles are positively 
charged (Z A = +210c, Z B = +100e) while after 48 hours 
the charge on the DiOCis-labeled spheres is reversed 
[Za = +110e, Zb — — 40e) and the two colloidal species 
are oppositely charged. Using the measured charges and 
an inverse screening length of k _1 = 400 nm, as deter- 
mined by conductivity measurements, we estimate the 
repulsive potential between unlike spheres at t = as 
~ +5/csT, at contact. After 48 hours, charge inversion 
generates a cross-attraction which is of order ~ — 3ksT. 
These estimates highlight the weak and reversible nature 
of the electrostatic interactions in our system. 

Hard-sphere suspensions are fluid at low volume frac- 
tions but crystallize into a random hexagonal close- 
packed structure at volume fractions (j> > 0.494. Crys- 
tallization continues up to a volume fraction 4> g « 0.58 
where a glass transition suppresses further nucleation and 
growth. A binary mixture of hard spheres of similar 
sizes is expected to be little different. For the size ratio 
7 = 0.93 used here computer simulations ^| and den- 
sity functional calculations confirm that the two species 
remain mixed in both the crystal and fluid phases and 
that there is no phase separation. To explore the effect 
of the cross attraction on the phase behavior we there- 
fore prepared binary suspensions with equal volumes of 
A and B particles and volume fractions which ranged 
from <fi = 0.48, just below the freezing transition up to 
4> = 0.58, the expected glass transition. The packing 
of the binary system of spheres was studied in real space 
with two-color fluorescence confocal microscopy. The dif- 
ferent emission spectra of DilCis and DiOCis dyes enable 
the two types of particles to be readily discriminated. A 
stack of 200 images, spaced by 0.16 fj,m vertically, was 
collected from a 73 /im by 73 /im by 32 /jm volume con- 
taining ~ 20000 particles of each species. The suspension 
was sealed in a cylindrical cell of ~ 50 /iL volume con- 
taining a short ~ 5 mm length of ferromagnetic wire to 
catalyze the decomposition of CHPB and which could be 
agitated with a magnet to rapidly mix the sample prior 
to observation. 

Suspensions with total volume fractions <fi < 0.487 



formed amorphous fluids with strong short-range order 
but no long-range correlations. An increase in volume 
fraction however led to the rapid formation of binary 
colloidal crystals. We found three distinct regular struc- 
tures depending on the degree of charging, as character- 
ized by the contact time t. Weakly-charged particles at 
high volume fractions assembled into a twelve-fold co- 
ordinated substitutionally disordered random hexagonal 
close- packed crystal (rhcp) (Fig. ^| . Computer simula- 
tions predict this to be a stable phase in hard sphere 
mixtures of similar size [l8| . In this structure the spheres 
are located at the crystal lattice sites in more or less ran- 
dom fashion. The lack of strong correlations between the 
two species is clearly evident in Fig. ^ 




FIG. 1: Confocal image of substitutionally-disordered rhcp 
lattice formed, after 23 hours, in a suspension with <j> = 0.507. 
The image depicts the close-packed {111} plane. The larger 
spheres (component A) are shown in green and the smaller 
spheres (B) in red. The inset figure shows the {111} face of a 
fee unit cell of binary spheres with volume fraction 4> c = 0.55. 

Increasing the magnitude of the opposite charges on 
each particle resulted in a completely different arrange- 
ment of large and small spheres. Fig.^a) shows a strik- 
ing two-color confocal micrograph of an AB superlattice 
formed in a suspension of = 0.528 after 199 hours. 
Fig. Hfb) depicts a crystallographic model of the su- 
perlattice built from a cesium chloride unit cell (space 
group SG 221). Rotation and cleavage along the {100} 
plane reproduced the structures observed in the confo- 
cal images. The cesium chloride structure consists of 
an equal number of cesium and chloride ions arranged 
alternately at the vertices of a body-centered cubic lat- 
tice so each particle has eight unlike neighbors. The lat- 
tice constant acsCi — 1.741 ±0.003 /^m, determined from 
Fig. |21 equates to a crystal volume fraction of <pcsC\ = 





FIG. 2: Confocal image and model of cesium chloride super- 
lattice (SG 221) formed, after 199 hours, in suspension with 
4> — 0.528. (a) Confocal image of {100} plane in polycrys- 
talline sample. Colors as in Fig. Note the square lattice 
of A spheres around each B. Inset shows CsCl unit cell, (b) 
and (c) Depiction of a slice through a CsCl crystal, parallel 
to the {100} plane, at three different depths of focus (dashed 
lines), (b) Projection perpendicular to {100} plane, (c) In- 
plane projection of {100} plane. Note the in-plane projection 
(c) remains unchanged as the depth of focus is stepped down 
through the crystal. This agrees with the experimental ob- 
servation that the symmetry of the confocal image (a) does 
not change as the depth of focus was varied. 



0.67 ± 0.01. The maximum packing limit for a CsCl 
crystal formed from a mixture of similarly-sized spheres 
( 7 > V3-1) is <^ sC1 (7) _= V3tt(1 + 7 3 )/[2(1 + 7) 3 ] ; For 
7 = 0.93 the corresponding maximum volume fraction is 
^CsCi = 0-68, which indicates that the CsCl superlattice 
observed is dense packed and the particles are close to 
touching each other. Note that although dense, the CsCl 
crystal is not a maximally-packed structure. When the 
two particles A and B are identically-sized, CsCl is iso- 
structural with a bcc packing of spheres which fills space 
significantly less efficiently than fee packing {<f)f cc = 0.74). 
Cesium chloride crystals were identified in binary samples 
with a wide range of volume fractions 0.507 < <fi < 0.578. 
In each case CsCl growth always proceeded initial for- 
mation of a substitutionally-disordered random-stacked 
close-packed crystal. 

A further increase in the charge on the two parti- 
cles led to the formation of a third AB superlattice. 
Fig. |3Ja) shows a confocal micrograph of the superlat- 
tice formed 215 hours after mixing in a binary mixture of 
total volume fraction <\> = 0.528. Crystallographic mod- 
eling (Fig. Gib)) reveals that the ordered regions visible 
are consistent with cleavage along the {100} and {111} 
(data not shown) planes of a NaCl superlattice (SG 225). 
The projection of the {100} plane visible in Fig. EJa) is 
superficially similar to the {100} plane of the CsCl super- 
lattice visible in Fig. |2Ia). However the two structures 
are easily discriminated by stepping the depth of focus. 



FIG. 3: Confocal image and model of sodium chloride super- 
lattice (SG 225) formed, after 215 hours, in suspension with 
4> = 0.528. (a) Confocal image of tilted {100} plane. Colors 
as in Fig.Q Note the square lattice of A spheres around each 
B. Inset shows NaCl unit cell, (b) and (c) Depiction of a slice 
through a NaCl crystal, parallel to the {100} plane, at three 
different depths of focus (dashed lines), (b) Projection per- 
pendicular to {100} plane, (c) In-plane projection of {100} 
plane. Projection (b) reveals that the {100} plane consists of 
vertical columns of alternate A and B spheres. Focusing down 
through the {100} plane (c) causes a shift in the 2D square- 
lattice of the face (shown in black) and a corresponding kink 
in the line of spheres. The same characteristic pattern is ev- 
ident in the confocal image (white dashed line) confirming 
that the square lattice seen in (a) arises from a NaCl rather 
than a CsCl structure. 



In the NaCl superlattice the vertically-aligned columns 
of particles, visible in Fig. [3J3, are composed of alternat- 
ing A and B particles in contrast to the CsCl structure 
(Fig. 12b). The NaCl structure is most conveniently envis- 
aged as two interpenetrating simple cubic lattices, occu- 
pied by an equal number of either sodium or chloride ions 
so that each ion has six unlike ions as its nearest neigh- 
bors. Analysis gave an average unit cell length aNaCi of 
3.3 ± 0.2 fim which equates to a crystal volume fraction 
of </>NaCi = 0.49±0.09. The maximum packing fraction of 
the NaCl superlattice is 0NaCi(7) = 2tt(1+7 3 )/[3(1+7) 3 ] 
for a mixture of spheres with a size ratio 7 > \/2 — 1. 
For 7 = 0.93 this expression gives a maximum volume 
fraction of </>NaCi = 0.53, which is in reasonable agree- 
ment with the value measured. The low packing fraction 
is readily understood if it is remembered that a NaCl 
crystal is equivalent to sc packing when no distinction is 
made between oppositely-charged spheres. It is striking 
that NaCl is observed even though 0NaCi ^ s considerably 
lower than the maximum packing fractions of either CsCl 
(fe cl = 0.68) or rhep (0 r * hcp = 0.74). 

The complex self-organization evident in our exper- 
iments is a consequence of a subtle interplay between 
entropy which favors dense close-packed structures and 
electrostatic interactions which prefer more open non- 
close packed structures with reduced particle coordina- 
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tion. The same balance plays a fundamental role in the 
study of molten salts where it is frequently described in 
terms of the restricted primitive model (RPM). The RPM 
consists of a mixture of equal-sized hard spheres, half of 
which carry a positive charge and the remainder an equal 
but opposite negative charge to ensure charge neutrality. 
The ions i and j interact through a potential which is a 
sum of a hard sphere repulsion and a purely Coulombic 
term. 



(r) 



oo 



if r < a 
f?1 if r > a. 



(1) 



In Eq. ^ Uc is the magnitude of the Coulombic potential 
at contact, a is the hard sphere diameter, and the plus 
and minus signs apply, respectively, to interactions be- 
tween ions of the same and different charge. The RPM 
is unlikely to be an accurate model of our experiments 
as it does not incorporate electrostatic screening. Never- 
theless the RPM accounts at least qualitatively for many 
of our observations. 

The solid phases formed at high density in the re- 
stricted primitive model have been inve stig ated by a 
number of authors, both theoretically Il9l 2fJ and more 
recently by Monte Carlo simulations |2ll 221 The 
phase diagram is a function of <f) and u c /fcsT, the rela- 
tive strength of the Coulombic potential at contact and 
the thermal energy. At low potentials u c < 3.5ksT, 
where hard sphere repulsions dominate, a fluid phase of 
ions is predicted to freeze at a density (j> w 0.5 into a 
substitutionally-disordered close-packed fee crystal. The 
transition is driven by the higher entropy of the fee crys- 
tal as compared to a fluid of the same density. Introduc- 
ing Coulombic interactions adds potential as well as en- 
tropic contributions to the total free energy. In a purely 
random fee structure each ion is surrounded by six ions 
of the opposite charge while in CsCl each ion has eight 
surrounding counterions. Consequently we expect a fec- 
CsCl transition as the Coulombic forces are increased. 
Computer simulations shows that the CsCl crystal is sta- 
bilized when u c > 4.3fcsT. At high densities and charges 
the RPM predicts a further solid phase [U . How- 
ever the structure is expected to be tetragonal rather 
than the cubic NaCl crystal observed here. Nonetheless 
it is well known that, for a number of alkali halides, the 
free energies of the CsCl and NaCl structures are suf- 
ficiently close that pressures of a few hundred kilobar 
are sufficient to inter-convert the two structures. Conse- 
quently our observation of NaCl may be a consequence 
of subtle differences in the interaction potential from the 
purely Coulombic mixture of oppositely-charged spheres 
assumed in the RPM. 

In summary we have demonstrated how the intro- 
duction of a weak reversible attraction between equal- 
sized particles leads to a surprisingly complex solid phase 
behavior with observation of substitutionally-disordered 
rhep, CsCl and NaCl binary crystals. Although to date 



we have explored only a single radius ratio (7 = 0.93) we 
have shown that a cross-attraction between particles sta- 
bilizes several new non-close packed structures, which are 
unstable in pure hard-sphere mixtures at the same size 
ratio. Our study suggests that competition between en- 
tropy and electrostatic interactions may provide a general 
route for the synthesis of new and more complex binary 
structures than can be stabilized by entropy alone and 
highlights the potential to direct the spontaneous growth 
of nano-particles into a diverse range of superstructures 
by engineering the interactions between species. 
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